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PROCEDE D' ESTIMATION D'UN RAYONNEMENT DIFFUSE , 
NOTAMMENT AFIN DE CORRIGER DES MESURES EN RADIOGRAPHIE 

DESCRIPTION 

5 Le sujet de cette invention est un procede 

d ' estimation d'un rayonnement diffuse, dont 
1' application principalement envisagee est la 
correction de radiographies . 

L 7 utilisation d'un rayonnement conique, 

10 tres fr^quente en radiographic, presente 1 ' inconvenient 
de produire un rayonnement diffuse important a travers 
l'objet examine. En d'autres termes, chacun des 
detecteurs situes derriere l'objet regoit non seulement 
un rayonnement primaire, proven'ant directement de la 

15 source par un trajet rectiligne et ayant traverse une 
region bien definie de l'objet, mais un rayonnement 
diffuse de provenance - indeterminee qui af f ecte la 
mesure et qu'il serait done souhaitable de corriger . 

Plusieurs procedes sont d6j^ pratiques. 

20 C'est ainsi que le rayonnement primaire peut §tre 
mesure seul si une collimation stricte des detecteurs 
et de la source est faite afin d' intercepter le 
rayonnement diffuse, mais ce procede necessite en 
pratique un balayage du faisceau qui est lent a 

25 accomplir, et pendant lequel on doit s'accommoder de 
mouvements du patient si on examine des etres vivants . 

On a aussi eu l'idee contraire de ne 
mesurer que le rayonnement diffuse. On dispose pour 
cela un reseau discontinu d' absorbeurs, coitune des 

30 billes de plomb, entre l'objet et les detecteurs, pour 
arreter local ement le rayonnement primaire, de sorte 
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que les detecteurs situes derri^re ces absorbeurs ne 
mesurent que. le rayonnement diffuse. Ce procede appele 
« beam stop » donne done des tables ou nappes 
bidiinensionnelles de valeur de rayonnement diffuse, 
5 qu'on complete par interpolation entre les detecteurs 
places derriere les absorbeurs. Le rayonnement diffuse 
ainsi estime est soustrait du rayonnement total mesure 
separement . Ce procede est precis mais a 1 ' inconvenient 
qu'il impose deux irradiations de l'objet et done un 

10 doublement de la dose de rayons qu'il regoit. Un 
dernier exemple de methode de correction du rayonnement 
diffuse par des moyens materiels comporte l'emploi de 
grilles anti-dif fusantes, mais leur efficacite n'est 
que partielle ; elle est insuf f isante pour un faisceau 

15 conique, ou le rayonnement diffuse peut etre plusieurs 
fois super ieur au rayonnement primaire . 

Enfin, il existe un certain nornbre de 
methodes numeriques pour estimer le rayonnement 
diffuse, a partir de convolutions ou de deconvolutions 

20 des mesures par exemple ; on pourrait aussi citer le 
brevet frangais 2 759 800 pour un procede numerique 
different, analytique* Elles sont en general d'emploi 
delicat car elles dependent de parametres choisis par 
1 'utilisateur (noyaux de convolution par exemple) qui 

25 ne donnent de bons resultats que dans des situations 
favorables, comme des petites zones ou le rayonnement 
diffuse est faible, ou des objets au contenu 
relativement homogene . Il n' existe aucun procede simple 
qui permette par exemple de corriger le rayonnement 

3 0 diffuse a travers le thorax ou d'autres grandes zones 
anatomiques , dont l'examen est frequent mais qui sont 
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defavorables pour corriger le rayonnement diffuse en 
raison de leur volume meme et de 1 'heterogeneite due a 
la presence d'une structure d'os complexe et ddnt la 
capacite d' attenuation du rayonnement est tres 
5 diff^rente de celle des tissus raous. 

Mentionnons enfin le brevet americain 
6 018 565 pour 1' expose d'une methode mixte, a « beam 
stop » et convolution. 

Un objet essentiel de 1' invention est de 
10 proposer un procede d' estimation et de correction de 
rayonnement diffuse qui puisse convenir pour des 
situations difficiles de radiographie . 

Le procede conforme a 1' invention est, sous 
sa forme la plus g^nerale, un procede d' estimation d'un 
15 rayonnement diffuse provenant d'un rayonnement initial 
ayant traverse un objet en subissant une attenuation 
laissant passer un rayonnement total de mesure, 
caracterise par : 

- une prise d'une table de mesures d'un 
20 rayonnement diffuse, obtenue en faisant passer le 

rayonnement initial par un simulacre de 1' objet, 

- un calcul de coefficients de 
transposition entre le simulacre et 1' objet, d'apres le 
rayonnement initial, le rayonnement total de mesure a 

25 travers 1' objet et un rayonnement total de mesure a 
travers le simulacre; 

- et une ponderation de la table de mesures 
avec les coefficients de transposition. 

Avantageusement , le simulacre sera un bloc 
30 d'epaisseur constante et en une matiere homogene, ayant 
une attenuation semblable a une matiere de base de 
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l'objet ; en general la prise de table de mesure sera 
une selection dans une s6rie de tables de mesures de 
rayonnement diffuse, obtenues auparavant en faisant 
success ivement passer le rayonnement initial a travers 
5 une serie respective de simulacres de l'objet, 
d'epaisseurs differentes mais constante ; et la 
selection sera faite par comparaison d'une valeur du 
rayonnement total de mesure a travers l'objet et d'une 
valeur du rayonnement total de mesure a travers les 
10 simulacres. 

Les coefficients de ponderation sont 
generalement des rapports de valeurs d'une meme 
f onctionnelle calculee pour l'objet et pour le 
simulacre. La f onctionnelle utilisee peut etre egale au 

15 produit du rayonnement total de mesure par le 
logarithme du rapport de rayonnement total de mesure et 
du rayonnement initial. 

L' invention sera maintenant. decrite en 
reference aux figures, parmi lesquelles : 

20 - la figure 1 est une vue generale d'une 

acquisition des mesures ; 

- la figure 2 est une vue d'une acquisition 
de calibration ; 

- et la figure 3 illustre les etapes du 

25 procede . 

Repor tons -nous d'abord a la figure 1, ou un 
tube 1 de rayons X emet un faisceau 2 conique vers un 
objet 3 a examiner (ici un patient etendu sur une table 
4) puis, a travers lui, vers un reseau 5 plan de 
3 0 detecteurs 6 disposes en matrice. Les detecteurs 6 sont 
relies a un appareil d' acquisition 7 et mesurent un 
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rayonnement diffuse qui se superpose au rayonnement 
primaire, seul convenable pour la radiographic. 

L' estimation du rayonnement diffuse a 
travers le patient 3 consiste tout d'abord a obtenir 
des tables bidimensionnelles ou nappes de rayonnement 
diffuse obtenues dans des circonstances comparables . 
Pour cela, on effectue des irradiations d'etalonnage a 
travers des simulacres 8 de l'objet 30 a examiner, 
conformement a la figure 2 : les conditions 
d' irradiation restent les memes, c'est-a-dire qu'on 
continue d'utiliser le tube 1, le faisceau 2, le reseau 
5 de detecteurs 6 et I'appareil d' acquisition 7, le 
simulacre 8 remplagant cependant le patient ; on a 
aussi ajoute une grille 9 de billes 10 de plomb entre 
le simulacre 8 et le reseau 5. II resulte de cette 
disposition que les rayons 11 passant par les billes 10 
sont completement absorbes et que les regions 12 du 
reseau 5 situees dans le prolongement de ces rayons 11 
ont des detecteurs 6 qui ne mesurent que le rayonnement 
diffuse a ces endroits. II suffit de relever ces 
valeurs mesurees et d'interpoler entre les regions 12 
pour es timer convenablement le rayonnement diffuse issu 
du simulacre 8 pour tous les detecteurs 6 du reseau 5 . 

Le simulacre 8 devrait etre semblable a 
l'objet afin que les rayonnements diffuses par eux 
fussent identiques. Une similitude parfaite h'est pas 
realisable, et c'est pourquoi on se contente d'un 
simulacre 8 ressemblant a l'objet 3 et dont la nappe 
associee de rayonnement diffuse sera corrigee 
ulterieurement pour evaluer celle de l'objet. En 
pratique, le simulacre 8 peut §tre un bloc d'une 
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matiere homogene et qui pr^sente le meme coefficient 
d' attenuation que la matiere de base de I'objet 3 : 
dans le cas d'un corps humain, compose pour l'essentiel 
de tissu mou, on sait que le plexiglas 
5 (polymethacrylate) convient. 

Afin de permettre des mesures variees , on 
disposera en realite de plusieurs nappes de rayonnement 
diffuse, obtenues pour autant de simulacres 8, qui ne 
differeront que par leur epaisseur et done par la 

10 longueur du trajet parcouru par les rayons 11, Ces 
nappes seront enregistrees dans une base de donnees 
prealablement aux mesures utiles sur les objets 3 
radiographics . pour prendre une nappe de rayonnement 
" diffuse comparable a celle d'un objet 3, on 

15 s^lectionnera en pratique une des nappes de la base de 
donnees ou, mieux, une nappe qu'on aura obtenue par des 
calculs d' interpolation entre deux de ces nappes. Le 
critere de selection pourra §tre defini au moyen d'un 
rayon particulier 13 aboutissant a une region 14 du 

20 reseau 5 et qui ne passera ni par les absorbeurs 10 de 
la figure 2, ni par des tissus osseux du patient (ou 
plus generalement des portions de I'objet 3 dont les 
proprietes d' absorption sont differentes du materiau du 
simulacre 8) a la figure 1. Le rayonnement total, 

25 primaire et diffuse, regu par la region 14 apres avoir 
traverse chaque simulacre 8 servira d' index a la table 
de rayonnement diffuse correspondante, et la table 
selectionnee aura 1' index a une valeur identique au 
rayonnement total mesurC a la region 14 a travers 

3 0 I'objet 3. Tout cela correspond au passage de l'etat El 
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a l'etat E2 dans 1 'organigramme de la figure 3, qu'on 
commence a coiranenter. 

La suite du procede consiste 

essentiellement en la correction de la table du 
5 rayonnement dif fusee ainsi selectionne pour 1'ajuster 
au mieux qu'on puisse esperer a la nappe de rayonnement 
reellement diffuse par l'objet 3. Pour cela, on se sert 
de toutes les informations disponibles, c'est-a-dire du 
rayonnement total regu par les detecteurs 6 au-dela de 
10 l'objet 3 comme du simulacre 8 selectionne. Ce 
rayonnement total etant note Ot, le rayonnement diffuse 
&d, le rayonnement initial issu du tube 1 <I>o et le 
rayonnement primaire <l>, la relation 4>t = <E> + ■ <I>d est 
respectee . 

15 On est alors aux etats E3 et E4 de 

1' organigramme de la figure 3. Ensuite, on transforme 
les valeurs des rayonnements totaux 4>t mesurees pour 
l'objet 3 et le simulacre 8 selectionne en leur 
appliquant des f onctionnelles . Plus precisement, il est 

2 0 connu dans 1'art que 4>d est proportionnel h 
4>log(4>/4>o) ; cette relation, qui est appelee la loi de 
Klein et Nishina, donne une allure generale du 
rayonnement diffuse, a defaut de son intensite. 

Le rayonnement initial <E>0 est connu ; le 

25 rayonnement primaire 4> ne I'est pas, mais on consent a 
appliquer cette relation de fagon approchee en le 
remplagant par le rayonnement total <£t, c'est-a-dire 
que la f onctionnelle employee associe a chaque valeur 
mesuree du rayonnement total Ot la valeur calculee 

30 Otlog(Ot/Oo) , supposee proche du rayonnement diffuse 
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Od k cet endroit ; on est parvenu aux £tats E5 et E6 de 
1 ' organigramme - 

L'etape suivante consiste a faire, pour 
chacun des detecteurs 6, le rapport des valeurs donnees 
5 par la f onctionnelle pour l'objet 3 et le simulacre 8 
selectionne selon la f ormule 

K = <^^^^blet Les coef f icients de 

Ot log(Ot / $o>imulacre 

ponderation K ainsi obtenus serviront a deformer la 
nappe de rayonnement diffuse selectionnee a l'etat E2 

10 afin d'es timer celle de l'objet 3. Les resultats 
constituent encore une table bidirnensionnelle ou une 
matrice ayant des dimensions identigues a celle des 
tables de rayonnement puisqu'elle est associee au 
reseau 5 de detecteurs 6* II est done possible et 

15 avantageux d'effectuer un filtrage numerique spatial de 
cette matrice en appliquant un filtre passe-bas qui 
corrige les coefficients K en ne conservant que les 
frequences les plus basses de leur variation et de les 
rendre - ainsi probablement plus conformes a la realite 

20 puisque le rayonnement diffuse varie assez lentement 
d'un point a un autre . 

Quand la table des coefficients de 
ponderation definitifs, notes K' , a ete obtenue (a 
l'etat E7), elle sert a ponderer la table de 

25 rayonnement diffuse selectionnee auparavant k l'etat 
E2, pour obtenir une table de rayonnement diffuse par 
l'objet 3 (etat E8 , qui constitue 1' estimation 
recherchee) ; la formule appliquee est 

Od objet=K' Od simulacre. Ces valeurs estimees 3>d 

3 0 objet pourront alors §tre soustraites du rayonnement 
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total 4>t mesur6 par les detecteurs 6 pour estimer le 
rayonnement primaire <J> et obtenir une image 
radiographique plus precise de l'objet 3. 

Ce procede s' applique aux radiographics a 
5 energie d' irradiation simple ou multiple ; dans le 
second cas, il est repute separ6ment pour chacune des 
energies employees . 

La f onctionnelle proposee ici n'est pas la 
seule qu'on puisse employer, et la f onctionnelle plus 
10 simple <&d=k<& (approchee ici encore en <t>d=k<I>t) , k etant 
une constante, pourrait aussi donner de bons r^sultats 
pour estimer <I>d. 

15 



B 13768*3 JCI 




diffuse provenant d'un rayonnement initial ayant 
traverse un objet (3) en subissant une attenuation 
5 laissant passer un rayonnement total de mesure, 
caracterise par : 

- une prise d'une table de rnesures d'un 
rayonnement diffuse, obtenue en faisant passer le 
rayonnement initial par un simulacre (8) de 1' objet, 

10 - un calcul de coefficients (K') de 

transposition entre le simulacre et 1' objet, d'apr^s le 
rayonnement initial ($o) , le rayonnement total de 
mesure a travers 1' objet (Ot objet) et un rayonnement 
total de mesure ci travers le simulacre (4>t simulacre) , 

15 - et une ponderation de la table de rnesures 

avec les coefficients de transposition. 

2. Procede d' estimation d'un rayonnement 
diffuse selon la revendication 1, caracterise en ce que 
le simulacre (8) est un bloc d'epaisseur constante et 

20 en une matiere homogene, ayant une attenuation 
semblable a une matiere de base de 1' objet. 

3. Procede d' estimation d'un rayonnement 
diffuse selon la revendication 1, caracterise en ce que 
la prise de table de rnesures est une selection dans une 

25 serie de tables de rnesures de rayonnement diffuse, 
obtenues en faisant successivement passer le 
rayonnement initial a travers une serie respective de 
simulacres de 1' objet, qui sont des blocs d'epaisseurs 
diff^rentes mais constante et en une matiere homogene, 

3 0 ayant une attenuation semblable a une matiere de base 
de 1' objet. 
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4. Procede d' estimation d'un rayonnement 
diffuse selon la revendication 3, caracterise en ce que 
la selection comprend une interpolation. entre deux des 
tables de mesures . 
5 5. Proced£ d' estimation d'un rayonnement 

diffuse selon I'une quelconque des revendications 3 et 
4, caracterise en ce que la selection est faite par 
coitiparaison d'une valeur du rayonnement total de mesure 
cl travers 1'objet et d'une valeur du rayonnement total 
10 de mesure a travers les simulacres , 

6. Procede d' estimation d'un rayonnement 
diffuse selon la revendication 5, caracterise en ce que 
la comparaison est faite pour des rayons identiques 
(13) du rayonnement initial a travers 1'objet et les 

15 simulacres, ne traversant que la matiere de base de 
1'objet. 

7. Procede d' estimation d'un rayonnement 
diffuse selon I'une quelconque des revendications 1 a 
6, caracterise en ce que les coefficients de 

2 0 ponderation sont des rapports de f onctionnelle 

identiques calculees pour 1'objet et pour le simulacre. 

8. Procede d' estimation d'un rayonnement 
diffuse selon la revendication 7, caracterise en ce que 
les f onctionnelles sont egales au produit du 

25 rayonnement total de mesure par le logarithme du 
rapport du rayonnement total de mesure et du 
rayonnement initial . 

9. Proced^ d' estimation d'un rayonnement 
diffuse selon I'une quelconque des revendications 1 a 

3 0 8, caracterise en ce qu'il comprend une etape de 

filtrage passe-bas des coefficients de transposition, 
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arranges en une table superposable a la table de 
mesures . 

10. Precede de radiographie comprenant une 
etape de correction de mesures de radiographie par une 
soustraction d'un rayonnement diffuse estime selon le 
procede selon 1'une quelconque des revendications 
prec^dentes . 
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